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Abstract 

Starting from the new background field formalism for T > 0, with 
nonperturbative (NP) background given as gauge invariant fiel corre- 
lators, new perturbation theory and diagram technic is introduced. 

Confined and deconfined phases are explicitly described and crit- 
ical temperature T c is expressed in terms of the scale anomaly term. 
Resulting numerical estimates of T c agree well with lattice data. Spa- 
tial area law is shown to follow naturally for T > T c , and a set of 
Hamiltonians is obtained for screening masses of mesons and glue- 
balla depending on the Matsubara frequences. 

Thelowest screening masses and wave functions of mesons obtained 
agree with lattice data and earlier calculations. Screening glueball 
masses and wave functions are also computed from the Hamiltonian 
and two different regimes are observed in high and low T regions. 

The infrared catastrophe of hot QCD is shown to be cured by the 
NP background. 



1 Introduction 

The finite temperature QCD provides a unique insight into the structure of 
the QCD vacuum, in particular it gives an important information about the 
mechanism of confinement. It is also of more practical interest, since the 
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nature of the deconflnement transition has its bearing on the cosmology and 
the hot QCD plasma can possibly be tested in heavy ion collisions. 

From theoretical point of view the hot QCD is a unique theoretical labo- 
ratory where both perturbative (P) and nonperturbative (NP) methods can 
be applied in different temperature regimes. 

It is believed that the perturbative QCD is applicable in the deconfmed 
phase at large enough temperatures T, where the effective coupling constant 
g{T) is small [1], while at small T (in the confined phase) the NP effects 
instead are most important. However even at large T the physics is not that 
simple: some effects, like screening (electric gluon mass), need a resummation 
of the perturbative series [2], while the effects connected with the magnetic 
gluon mass demonstrate the infrared divergence of the series [3]. 

During last years there appeared a lot of lattice data which point at the 
NP character of dynamics above T c . Here belong i) area law of spacial Wilson 
loops [4] ii) screening masses of mesons and baryons [5] and glueballs [6] iii) 
temperature dependence of Polyakov-line correlators [7] . 

In addition, behaviour of e — 3p above T c has a bump incompatible with 
the simple quark-gluon gas picture [8]. 

Thus the inclusion of NP configurations into QCD at T > and also at 
T > T c is necessary. 

Recently a systematic method for QCD was developed, treating NP fields 
as a background and doing perturbative expansion around that both for 
T = [9] and T > [10]. 

To describe the phase transition a simple choice of deconfmed phase was 
suggested where all NP color magnetic configurations are kept intact as in the 
confined phase, whereas colorelectric correlators responsible for confinement, 
vanish. 

This picture together with the background perturbation theory forms a 
basis of quantitative calculations, where field correlators (condensates) are 
used as the NP input. 

The plan of the lecture is as follows. In the second chapter the new 
background field formalism is presented, based on the familiar background 
field method [11] modified for T > and incorporating the 'tHooft's identify 
for integration over quantal and background fields. 

The temperature phase transition is discussed in chapter 3 and resulting 
predictions for T c are compared with lattice data. 

In chapter 4 the spacial Wilson loops are computed both below and above 
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T c ; it is shown how at large T a new regime - that of dimensional reduction 
- appears and the spacial string tension is discussed. 

Chapter 5 is devoted to screening masses of mesons, and glueballs, which 
were considered recently in [5,6]. 

In chapter 6 the analysis is given of the infrared catastrophe of the hot 
QCD; it is shown, that it is cured naturally in the present formalism by 
nonperturbative contributions. 

Some prospectives of the method and discussion of results are given in 
Conclusions. 



2 New background field formalism 

We derive here basic formulas for the partition function, free energy and 
Green's function in the NP background formalism at T > [10]. The total 
gluonic field is split into a perturbative part and NP background 

An = B„ + a„ (1) 

where both and subject to periodic boundary conditions. The 

principle of this separation is immaterial for our purposes here, and one can 
average over fields B^ and a M independently using the 'tHooft's identity^] 

Z = / DA^ P (S(A)) = IDB^BU^-SV + a)) p) 

=<< exp(—S(B + a) > a >B 

with arbitrary weight T)(B). In our case we choose rj(B) to fix field correlators 
and string tension at their observed values. 
The partition function can be written as 

Z(V,T,n = 0) =< Z(B) > B , 
Z(B) = N J D<f)exp(- dr J d 3 xL(x, r)) (3) 



1 private communication to one of the authors (Yu.S.), December 1993. 
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where denotes all set of fields a M , N is a normalization constant, 

and the sign <>b means some averaging over (nonperturbative) background 
fields B^, as in (2) . Furthermore, we have 

8 

L(x,t) = ^2 Li, 

i=l 

where 

L 1 = \(F^(B)) 2 ;L 2 = ^W^, 
L 3 = Q a (D 2 (B)) ab Q b ; L 4 = -igQ a (D^, a^ ab Q b (4) 
L 5 = -^a(D ll (B)a ll ) 2 ;L 6 = L int (a 3 ,a*) 
L 7 = —a v D tl (B)F flv (B); L 8 = ^+(m + D(B + a))* 

Here 6, 6 are ghost fields, a- gauge-fixing constant, L 6 contains 3-gluon- 
and 4-gluon vertices, and we keep the most general background field B^, not 
satisfying classical equations, hence the appearance of L 7 . 

The inverse gluon propagator in the background gauge is 

W; b u = -D\B) ab • V - 2gF; v {B)r b (5) 

where 

(D x ) ca = d x 5 ca - %gT h ca B\ = d x 5 ca - gf bca B b x (6) 

We consider first the case of pure gluodinamics, L 8 = 0. 
Integration over ghost and gluon degrees of freedom in (3) yields 

Z(B) = N\detW(B)); e f[det(-D,(B)D,(B + a))] a=i _ x 

oo ql x 1 

x {l + J2^f(a = J1 )}exp(--JW- 1 J) J ^ d^f^b) (7) 

One can consider strong background fields, so that gB^ is large (as com- 
pared to Aq CD ), while a s = in that strong background is small at all 
distances [9]. 

In this case Eq. (7) is a perturbative sum in powers of g n , arising from 
expansion in (ga^) n . 



4 



In what follows we shall discuss the Feynman graphs for the free energy 
F(T), connected to Z(B) via 

F(T) = -Tin < Z(B) > B (8) 

As will be seen, the lowest order graphs already contain a nontrivial 
dynamical mechanism for the deconfinement transition, and those will be 
considered in the next section. 

To the lowest order in ga^ the partition function and free energy are 

Z =< exp(-F (B)/T) > B , 
F (B)/T = hndetW - lndet{-D 2 (B)) = (9) 

= S P [ at)j{-\e- tW + e^) 

where W = —D 2 (B) — 2gF and D 2 (B) is the inverse gluon and ghost prop- 
agator respectively, ((t) is a regularizing factor [10]. 
The ghost propagator can be written as [10], [12] 

(-DX 1 =< x\ j™dte tD2 ^\y >= ^ dt(Dz): y e- K $(x,y) (10) 

where standard notations [12] are used 

and a winding path integral is introduced [10] 

(Dz) w - lim TT V f ^ MTx( m )-(x-y)-nps ui ) , n s 



with /3 = 1/T. For the gluon propagator an analogous expression holds 
true, expect that in (4) one should insert gluon spin factor PpexplgF inside 
y). For a quark propagator the sum over windings in (5) acquires the 
factor (— l) n and quark spin factor is expga^F^ [10]. 
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3 The temperature phase transition in QCD 



We are now in position to make expansion of Z and F in powers of ga^ (i.e. 
perturbative expansion in a s ), and the leading-nonperturbative term Zq,Fq 
- can be represented as a sum of contributions with different N c behaviour of 
which we systematically will keep the leading terms 0(N 2 ), 0(N C ) and 0(iV°). 

To describe the temperature phase transition one should specify phases 
and compute free energy. For the confining phase to lowest order in a s 
free energy is given by Eq.(3) plus contribution of energy density e at zero 
temperature 

F(l) = eV 3 - |V 3 T 4 -T^ ^ffi e^ + Q(l/N e ) (12) 
where e is defined by scale anomaly [13] 

£ ~ ~ Y Nc ^t < {F ^ B)f > (13) 

and the next terms in (12) correspond to the contribution of mesons (we keep 
only pion gas) and glueballs. Note that e = 0(N 2 ) while two other terms in 
(12) are 0(iV c °). 

For the second phase (to be the high temperature phase) we make an 
assumption that there all color magnetic field correlators are the same as in 
the first phase, while all color electric fields vanish. Since at T = color- 
magnetic correlators (CMC) and color-electric correlators (CEC) are equal 
due to the Euclidean 0(4) invariance, one has 

< (Fpu(B)) 2 > = < (F£ v ) 2 >el + < (F°„) 2 >magn] < F 2 > m agn = < F 2 > ei 

(14) 

The string tension a which characterizes confinement is due to the electric 
fields [14], e.g. in the plane (i4) 

2 

a = a E = i- J J d 2 x< trEi(x)$(x,0)Ei(0)<i>(0,x) > +... (15) 

where dots imply higher order terms in Ei. 

Vanishing of <Je liberates gluons and quarks, which will contribute to 
the free energy in the deconfined phase their closed loop terms (10) with all 
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possible windings. The CMC enter via perimeter contribution 

< $(x, x) >= Q (see (9,10)). As a result one has for the high-temperature 

phase (phase 2) (cf.[10]). 



F(2) = \eV, - (N 2 - l)V^Q g - 7 -^N c V 3 T%Q q + 0(N° C ) (16) 



where Q q and Q g are perimeter terms for quarks and gluons respectively, the 
latter was estimated in [15] from the adjoint Polyakov line; in what follows 
we replace Q by one for simplicity. 



Comparing (12) and (16), F(l) = F(2) at T = T c , one finds in the order 
0(iV c ), disregarding all meson and glueball contributions 



For standard value of G 2 = ^ < F 2 >= 0M2GeV A [13] (note that for 
rif = one should use approximately 3 times larger value of G 2 [13]) one has 
for SU (3) and different values of rif — 0, 2, 4 respectively T c = 240, 150, 134 
MeV. This should be compared with lattice data [8] T c (lattice) = 240, 146, 131 
MeV. Agreement is quite good. Note that at large N c one has T c = 0(N®) 
i.e. the resulting value of T c doesn't depend on N c in this limit. Hadron 
contributions to T c are 0(N~ 2 ) and therefore suppressed if T c is below the 
Hagedorn temperature as it typically happens in string theory estimates [16]. 

Till now we disregarded all perturbative and nonperturbative corrections 
to F{2) except for magnetic condensate, the term \sV,. If we disregard also 
this term, considering in this way only free gas of gluons and quarks for the 
phase 2, we come to the model, considered in [17]. The values of T c obtained 
in this way differ from ours not much - they are factor of 2 1 / 4 = 1.19 larger, 
but one immediately encounters problems with explanation of spacial string 
tension, screening masses etc., which are naturally are accounted for by the 
notion of magnetic confinement - nonzero values of magnetic correlators in 
the phase 2, including magnetic condensate term, ^sV 3 . 

Our approximation (10) corresponding to lowest order in Nc and g was 
too simplified when we have put fl g = fl q = 1. 

Indeed NP corrections may contribute to Q g , Q q . Their phenomenological 
necessity can be seen in the measured values of e — 3p, [7,8], which are seen 




(17) 
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in Fig.l In case of £l g = Q q = 1 the difference e — 3p should be zero, and of 
course higher orders in N^j 1 (NP effects ) and higher orders in g (Perturbative 
effects) contribute to it. In [15] we have tried to estimate effect of nonzero 
(Q g — 1), which is O(JV^), on the energy density and preassure. To this end 
we exploit the adjoint Polyakov line and separate from it the NP perimeter 
contribution. This can be done if one subtracts properly the linear divergent 
perturbative contribution, specific for Wilson and Polyakov contours. Thus 
one can write [15] 

^ (iVP) = X(pert ' n9iPert) = eX ^-f L mNc] (18) 

Substituting this value of Q g (NP) into F(2), eq.(16) , one obtains e and p, 
shown in Fig. 2 by the line. One can see a reasonable agreement between 
this estimate and lattice data, which could signify the importance of gluon 
perimeter contribution. That was done for the SU (2) group, since only there 
exist detailed lattice data for the adjoint Polyakov line. 

One can now estimate the influence of nonzero Q g — 1 on Tc, again for 
the SU (2) case. 

Taking Q g (NP) from (18) and substituting it into (17) one obtains a shift 
of T c by a factor of 1.07 for rif — 2 [18]. Similar calculations for SU(3) are 
now in progress [18]. 



4 Spacial Wilson loops 

In this section we derive area law for spacial Wilson loops, expressing spacial 
string tension in terms of CMC. 

To this end we write < W(C) > for any loop as [14] 

2 

< W(C) >= ex V \- 9 — J da^{u)da pX {u) < F^(u)<$>(u, u')F pX (u')<S>(u' , u) > 

(19) 

+ higher order cumulants] 

For temporal Wilson loops, in the plane i4, i = 1,2, 3, only color electric 
fields Ei = Fj 4 enter in (19), while for spacial ones in the plane i,k;i,k = 
1,2,3 there appear color magnetic field Bi = \eikiFki; in standard way [14] 
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one obtains the area law for large Wilson loops of size L, L ^> (Tj m ) is 
the magnetic correlation length) 

< W(C) > sp aciai~ exp(-a s S min ) (20) 
where the spacial string tension is [10,14] 
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J d 2 x < B n (x)$(x, 0)£ n (0)$(0, x) > +0(< B 4 >) (21) 



and n is the component normal to the plane of the contour, while the last 
term in (21) denotes contribution of the fourth and higher order cumulants. 
On general grounds one can write for the integrand in (21) 

dD B dD B 

< B i (x)^(x,0)B J (0)^(0,x) >= S ij (D B (x)+D B (x)+x 2 -^j-)-x i x j -^, 

(22) 

and only the term D B (x) enters in (21) 
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J d 2 xD B (x) + 0(< B A >) (23) 



similarly for the temporal Wilson loop in the plane i4 one has the area law 
for T <T C with temporal string tension 

2 

a E = Yf d 2 xD E (x) + 0(< E 4 >) (24) 

For T = due to the 0(4) invariance CEC and CMC coincide and oe = cr s . 
For T > T c in the phase (2) CEC vanish, while CMC change on the scale of 
the dilaton mass ~ lGeV, therefore one expects that a s stays intact till the 
onset of the dimensional reduction mechanism. This expectation is confirmed 
by the lattice simulation - a s stays constant up to T 1.4T c [19]. Recent 
lattice data [19] show an increase of a s at T ?a 2T C , for SU(2) which could 
imply the early onset of dimensional reduction. 

Indeed the numerical analysis [7] shows that at T > 2T C the spacial string 
tension a s can be well reproduced as 



^{T) = cg 2 (T)T (SU 3 ) (25) 
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which should be compared with the 3d QCD value [20] 



v^T = 0.55% 3 2 (ST/3) (26) 

The constant c in (25) is actually expandable as a series of g n , but in 
the range 2 < T/T c < 4 it is approximately constant, c ~ 0.63. Comparison 
of (25) and (26) with g\ = g 2 (T)T indeed supports an early dimensional 
reduction with the coupling constant in the 2-loop approximation 

-aL = Wn(-^) 2 + ^Inln^f (27) 
g 2 (T) A T b A T 

where 60 = 1^2, &i = (ilpp and the value At fitted to cr s (T) is 

A- = (0.076 ± 0.013)T C (28) 

Due to the small value of Ax,g(T) is also small, g(2T c ) as compared 

to the 4d, T = value from heavy quarkonia, A ps 200MCV. 

Question: Why dimensional reduction sets in as early as T = 2T C ? 



5 Screening masses of mesons and glueballs 

This section is based on results of ref. [21]. In this section we consider the qq 
and gg Green's functions G(x, y) at T > T c and derive corresponding Hamil- 
tonians for evolution in the spacial direction. We start with the Feynman- 
Schwinger representation [22] for G(x, y), where we neglect for simplicity spin 
interaction terms 

G(x,y) = J o ds J o dse- K - K {Dzy xy {DzY xy < W{C) > (29) 

Here K and (Dz)™ y are defined in (10) and < W{C) > in (20), where the 
contour C is formed by paths z(r), z(t) and t = x — y is for definiteness 
along the axis 3. Since by definition at T > T c electric correlators are zero, 
only elements da^ u in (19) in planes 12,13 and 23 contribute. As a result one 
obtains for < W(C) > the form (20) with 

Smin = dr d'y^Jwfw'^ - (w^) 2 (30) 
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where only spacial components i = 1, 2, 3 enter 



dzv ■ dw ■ 

Wi(r,7) = ^(r)7 + ^(r)(l -7),^ = , w • = (31) 

The form (30) is equivalent to that used before in [12] but with w 4 = 0. 

As a next step one can introduce " dynamical mass" fj,, p similarly to 
[12]. We are looking for the "cm" Hamiltonian which corresponds to the 
hyperplane where z 3 = z 3 . Now the role of evolution parameter (time) is 
played by z 3 = z 3 = r with < r < t, and we define transverse vectors 
z±(z 1 ,z 2 ),z±(z 1 ,z 2 ) and z 4 (t), z 4 (r). 

dz 3 dr dz 3 dr 

lx=Tx = 2 ^ -dx=Tx= 2 ^ (32) 

then K 7 K in (29) assume the form 

1 ft 771 ^ 

K = g I dr[-^ + Mr)(l + z\ + z\)\ (33) 

and the same for K with additional bars over /i, z±,z 4 . 

Performing the transformation in the functional integral (29) dsDz 3 {r) — > 
Dfi, dsDz 3 {r) — > Dp one has 

G(x,y) = J DfiDpDz ± Dz ± (Dz 4 )™ y (Dz 4 )™ y exp(-A) (34) 

with the action 

A = K + K + aS min (35) 

Note that z 4 , z 4 are not governed by NP dynamics and enter A only kinemat- 
ically (through K,K), and hence can be easily integrated out in (34) using 
Eq.(ll) for the 4-th components - with (x — y) 4 = 0. One can now proceed 
as it was done in [ 12], i.e. one introduces auxiliarly functions u(r, 7), t](t, 7); 
defines center-of-mass and relative coordinates R±, f± = f, and finally inte- 
grates out R± and t](t, 7). The only difference from [12] is that now z 4 , z 4 do 
not participate in all those transformations. As a result one obtains 

G(x,y) ~ f DvD^DpDre~ A(1) J2 e~ An ^ (36) 



n,n2 
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here 

^2 i ™2 ^2 i ™2 



A {1 > /i, /i, i/ = - / dr - + - _ 2 +fi + fL + a 2 r 2 / — + / vdi+ 

2 Jo /x /i Jo is Jo 

(37) 

L 2 /r 2 



Ml-C) 2 + K 2 + /o^7(7-C)V 

A — = 2 (7rT) y rfr( MT + (38) 

b(n) = 2n for bosons and 2n + 1 for quarks. We also have introduced radial 

— # 

momentum p r , angular momentum L 



fr s <EI = ( J^<r3! f = fX p (39) 



and 

„(r)+< ? >J^ < 7>s i^2 (40) 

Let us define the Hamiltonian for the given action 
in (36), integrating over Du, D/j,, Dp, around the extremum of A ( this is an 
exact procedure in the limit t — > 00). For the extremal values of auxiliary 
fields one has 

f + m \ + (fo^OTrT) 2 _ _ /(/ + !) (1-C) 2 _ 

1 ^> V „2 ..9! ^ 



/i 2 (r) f 2 a§ /i : 

^ + m 2 + (6( n2)vr T) 2 _ x _ 1(1 + 1) C 2 _ J_ } (41) 



/i 2 (r) a 2 /i 2 

r = 1 — 



^ 1 ^ + 1)(7"0 2 



z/ 2 (r, 7) r 2 a 2 

where 03 = — £) 2 + /X 2 + / ^7(7 ~~ C) 2// an d C is defined by eq.(40). 

After the substitution of these extremal values into the path integral 
Hamiltonian 

G(x,y) =< x\ e~ Hn ^ l \y > (42) 

one has to construct (performing proper Weil ordering ) the operator Hamil- 
tonian acting on the wave functions. 
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Consider for simplicity the case L = 0, then from (37-38) one obtains 



H ni n 2 = sjv 1 + m\ + (&(ni)7rT) 2 + ^ + m\ + {b^Tf + or (43) 

Here p = 4Jj and f — f± is a 2d vector, f— (ri, r 2 ); mi, m 2 - current masses 
of quark and antiquark, for gg system mi = m 2 = and a = a a dj = \c ■ 
Eigenvalues and eigenfunctions of H nin2 

H nina ip(r) = M( ni ,n 2 )^(r) (44) 

define the so-called screening masses and corresponding wave-functions, which 
have been measured in lattice calculations [5,6,23]. 

The lowest mass sector for mesons is given by H 00 , where for n\ = n 2 = 
one has b{n-y) = b(n 2 ) = 1 in (43). For light quarks one can put m x = m 2 = 
and expand at large T square roots in (43) to obtain 

#oo ~ 2vrT + L- + ar, M 00 ~ 2vrT + e(T) (45) 

where e(T) = Q and a ~ 1.74 is the eigenvalue of the 2-D dimension- 

less Schrodinger equation. 

Assuming the parametrization ^J~o~ s = cg 2 (T)T with c = 0.369 and scaling 
behaviour of g 2 (T) [19] one has M 00 ~ 2nT + 0((ln(T/A T ))" 4 / 3 T) tending 
to twice the lowest Matsubara frequency. (This limit corresponds to the 
free quarks, propagating perturbativly in the space-time with the imposed 
antiperiodic boundary conditions along the 4th axis). 

Eq. (45) coincides with that proposed in [24] , where also numerical study 
was done of M o and ipoo( r ) ■ Our calculations of Eq. (44-45) agree with 
[24] and are presented in Fig.l together with lattice calculations of tpoo(r) 
for p - meson [23]. The values of M o(T) found on the lattice [3,4] are 
compared with our results in Fig. 4. Note, that our M o (45) does not contain 
perimeter corrections which are significant. Therefore one has to add the 
meson constant to M o to compare with lattice data. We disregard spin- 
dependent and one-gluon-exchange (OGE) interaction here for lack of space. 
It is known [24] however, that OGE is not much important at around T ps 
2T C . 
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We note, that the lowest meson screening mass appear also in Yukawa 
type exchanges between quark lines and can be compared with the corre- 
sponding Yukawa potentials. 

For the gg system the lowest mass sector is given by b(n = 0) = 0, and 
one has from (43) 

H 00 (gg) = 2\p\ + a adj r (46) 

Note that T does not enter the kinetic terms of (46). 

To calculate with (46) one can use the approximation in (37) of r- 
independent fi [12], which leads to the operator (/i = fl) 

f 

h(fii) =n-\ h cr adj r (47) 

fx 

The eigenvalue E(fi) of h(fi) should be minimized with respect to fx and 
the result E(fx ) is known to yield eigenvalue of (46) within few percent 
accuracy [25]. The values E(fx ) xs M gg thus found are presented in Fig.4. 
The corresponding wave functions ipoo(r) are given in Fig.3. These data can 
be compared with the glueball screening masses, found on the lattice in [6]. 
Another point of comparison is Polyakov line correlator 

P(R) =< n(R)Q(0) > - < Q(R) >< ft(0) >= exp{-V{R)/T) - 1 (48) 

with 

V(R) = eXp{ '^ R) , Q(R) = -LtrPexptg f A 4 dz 4 
R a JNc -Jo 

It is easy to understand that fi = M gg and one can compare our results for 

glueball screening mass (GSM) M gg (T) with the corresponding lattice data 

for fx in Fig.4. 

Taking again into account an unknown constant perimeter correction to 
our values of M gg , one can see a reasonable qualitative behaviour. 

In addition to the purely nonperturbative source of GSM described above 
there is a competing mechanism - the perturbative formation of the electric 
Debye screening mass m e i(T) = gT{^- + ^) for each gluon with g(T) given 

by the temperature dependence of ^a s {T) = cg 2 (T)T [7,19]. Therefore for 
large T, where m e i is essential, one should rather use instead of (46) another 
Hamiltonian, which is obtained from (43) replacing mi = m 2 = m e i{T) 



H 00 (gg) = 2JfP + m 2 el (T) + a adj (T)r (49) 
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It turnes out that up to the temperatures T < 2T C one can consider the effect 
of m e i(T) as a small correction to the eigenvalue of (46) E 

M gg = E + 5E* ■ \cg\T)T + (50) 

and that gives us the right to treat GSM in the temperature region concerned 
nonperturbatively. 

We note that the transition between these two regimes (with the domi- 
nance of nonperturbative and then perturbative dynamics) is smooth, with 
GSM tending to twice the Debye mass at large T. 



6 Infrared catastrophe in Hot QCD 

The infrared (IR) catastrophe in QCD was indentified within the hot per- 
turbative series in [3]. For our purposes we shall explain it in the language 
of Feynman diagrams in the configuration space. Consider the Feynman di- 
agram for the thermodynamic potential in the order g n , shown in Fig. 5. It 
can be written as 

n 

jW = l[d 4 x^r^l[G(x^,x^) (51) 

i=l i,j 

where Gij is the x- space gluon Green's function connecting vertices % and 
j, and is the three-gluon vertex containing one derivative, d/dx®. We 
suppress all color and Lorentz indices for the sake of brevity and will only 
count the powers of to judge the IR convergence of J( n \ 
At T > using Matsubara formalism one has to replace 

d 4 x -> T^cPx 

/Trfin p —w(2—y)-i2TrkT(x4,—y4,) 
_ id ( P -+w> < 52 » 

For large T the only mode propagating over distances larger that 1/T, as 
follows from (52) is the mode with k = 0. Exactly this mode causes IR 
catastrophe; we have for it 

G(x,y)~ J |f-y|»l/T (53) 

47i\x — y\ 
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Now for the diagram of the type of Fig. (5) one has for a given number of ver- 
tices n total number of propagators ^p, so that one can write (remembering 
that one vertex yields overall volume) 



J K > ~ g V 3 1 

x 2 

~ g n V 3 T* J (xf- 3 ) 



(54) 



where the the symbol J(x k ) denotes the overall power of in the integrand 
including those of d 3 x^\ One can immediately see that for n > 6 the diagram 
diverges at large x. 

We shall now give simple arguments which show that accounting for the 
nonperturbative background, the perturbation series for F behaves well and 
IR catastrophe disappears. 

To this end we consider the perturbation theory given by Eqs. (1-11). 

The diagram of Fig. (5) goes over into the equivalent diagram of the order 
g n , where vertices and propagators are now taken in the NP background field. 

[ ' V) \ (Wy x l y = (-D 2 - 2gF)£, for the gluon [bb) 
rW ~ d/dxf - - igB^x®) 

OXfi 

The expression for the diagram, looks the same as (53) with eq.(10), for 
propagators, but now two major changes appear: 

i) All phase factors <f>(x, y) of gluon (ghost) Green's function and vertixes 
T assemble into gauge-invariant combinations of Wilson loops and covariant 
derivatives, so that the diagram is gauge-invariant. 

ii) Averaging over background fields as in (2) brings about Wilson-loop 
averages, which are subject to area law for large loops (as can be seen from 
cluster expansion [14] or lattice calculations [19]). 

Therefore it is highly plausible that large distance behaviour of the to- 
tal diagram is cut - off at large distances \x® — by the Wilson loops; 
by dimentional arguments it reduces to saying that the divergent integral 
f(x%~ 3 ), (n > 6) should be replaced by fx 3 ~ n ^ 2 , with ii m y/o^- 
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Correspondingly Eq.(54) is replaced in the framework of the NP back- 
ground perturbation theory by 

j(») ~ g n VzT n l\^F:f- n W (56) 

Now at high T, T 3> 2T C , as we discussed in chapter 4, the dimensional 
reduction mechanism is at work, which yields [7,19] ^a s {T) = cg 2 (T)T. 
Inserting this into (56) one gets 

jW rsj V 3 T 3 (cg 2 (T)) 3 (^-) n/2 (57) 

c 

Thus all diagrams with n > 6 are not IR diverging due to the Spatial 
string tension a s and their contribution to the free energy density can be 
written in the same form as in [26] 

^ = ^ = {9 2 Tff G + ... (58) 

where fa = J2n>6 a n(^) 11 -^ is given by the sum of diagrams with n > 6, 
and can be easily computed since g(T), Eq.(27) is small for large T. Note 
however, that c = c(g) as discussed, below Eq.(26). 

Since our conclusion about the absence of the IR catastrophe is crucial for 
the theory, we present now another argument on convergence of JH To this 
end let us consider behaviour of the diagram, Fig.5., when one of the vertices, 
say i — 1, is far from all others. The diagram is actually longdistance limit 
of 3g glueball Green's function, since 3 gluons are produced at the vertex 1 
in the gauge-invariant state with total angular momentum zero (this vertex 
is given by Lq(o, 3 ) term in (4)). 

It is clear, that at large \x^ — X\, where X stands for the assembly of 
all other vertices, the diagram behaves as 

jW(x«,x( 2 L.) ~ e -M 3g (o+)\ x m-x\ (5Q) 

where M 39 (0 + )is the mass of the 3-gluon glueball in 3d, one expects that 
M 3g (0+) > lGeV. 

Now M 3g (0 + ) ~ y/o^const, where const is some number of order of one. 
Thus again we obtain a cut-off in the x-space, which makes all integrals 
convergent and the estimate (57) remains true. 
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We note that nonperturbative background has created a situation equiv- 
alent (in the sense of convergence) to the appearance of the magnetic mass, 



Note however, that M magn appears in a gauge-invariant manner, and not 
in Yiij ~ self energy part of gluon. Rather it appears as in the constituent 
gluon model; if one defines M gg = 2M magn (cf.Eq.(50)), then M 3g 3M magn , 
and each gluon propagator contributes in x-space assymptotics G(x, y) ~ 
e -M magn \x~y\^ ^his asymptotics reproduces immediately (59) and (56). Also 
the Polyakov line correlator (48) containing the exchange of two gluon lines, 
decreases at large distances with screening mass \i = M gg ss 2M magn . 



We have presented the formalism of QCD at T > 0, which contains NP field in 
the background, in the form of gauge-invariant correlators and perturbative 
expansion in this background. 

It is shown how the Feynman-Schwinger representation (FSR) which 
proved to be very useful for T = 0, is modified for T > 0. In particular, 
all gluon and quark propagators are written as a sum of path integrals with 
winding paths in the 4th direction. 

Using this formalism the structure of vacuum at T < T c and T > T c as 
well as the nature of deconfinement transition is clarified and T c is computed 
in good agreement with lattice data. Moreover the spacial area law is shown 
to follow naturally from the structure of the vacuum at T > T c . Using FSR 
the Hamiltonian for screening states of mesons and glueballs is derived, and 
screening masses and wave functions for lowest meson and glueball states are 
numerically computed. 

Finally, it is demonstrated how the Infrared Catastrophe is eliminated 
due to the NP background and estimates of higher order diagrams for the 
free energy are given. All this exemplifies NP background field theory and 
in particular perturbation theory as a powerful instrument both for T <T C 
and for T > T c . 

The author acknowledges the financial support from the Russian Fund 
for Fundamental Research, Grant 93-02-14937. 
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Figure captions 

Fig. 1. The intraction measure, (e — 3p/T 4 as a function of temperature, 
from ref. [7]. 

Fig. 2. The normalized pressure p/psb and energy density e/esb calcu- 
lated from (16), (18) (p. back., e.back) and lattice calculations (p.latt., e.latt) 
as functions of T/A L , from ref. [15]. 

Fig. 3. The p- meson (solid line) and glueball (dashed line) wave, func- 
tions for lowest Matsubara frequencies vs r/a(a = 0.23/m) for T = 210MeV 
from ref. [21]. The lattice data are from ref [23] for the same T. 

Fig. 4. The screening masses for mesons (solid line) and glueballs (dashed 
line) as functions of temperature from ref. [21]. The lattice data are from 
ref. [5] (squares) and ref. [6] (triangle). 

Fig. 5. The typical IR divergent diagram of hot QCD in the order g n . 
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Fig.l, Pig.2, Pig.3, Fig.4, Fig.5, 



^3? ^4? ^n—l) %n 
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